Unified Brane Gravity: 
Cosmological Dark Matter from Scale Dependent Newton Constant 
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We analyze, within the framework of unified brane gravity, the weak-field perturbations caused 
by the presence of matter on a 3-brane. Although deviating from the Randall-Sundrum approach, 
the masslessness of the graviton is still preserved. In particular, the four-dimensional Newton force 
law is recovered, but serendipitously, the corresponding Newton constant is shown to be necessarily 
lower than the one which governs FRW cosmology. This has the potential to puzzle out cosmological 
dark matter. A subsequent conjecture concerning galactic dark matter follows. 
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1. Introduction 

The discovery of dark matter [l| continues to be a 
challenging problem for astrophysics and cosmology. Al- 
though many ideas from particle physics have been 
put forth, none of them so far have been able to provide 
a convincing explanation for its mysterious nature and its 
preponderance in the constitution of the Universe. Be- 
sides, the search for a suitable particle candidate so far 
has proved elusive. A different approach proposes that 
dark matter is not matter at all in the conventional sense 
but rather an artifact originating from a deviation from 
general relativity. Various forms of modified theories of 
gravity have been proposed to explain the phenomenon 
[3|. However, all such theories are beset with their own 
problems. In this context, it is important to note the re- 
markable coincidence in the amounts of galactic and cos- 
mological dark matter. This is natural for particle dark 
matter but has to be explained by any full theory of arti- 
fact dark matter in modified theories of gravity. Recently, 
the idea that brane world theories could provide an expla- 
nation for the dark matter has been suggested Q . Brane 
world theories have recently made great breakthroughs 
in the area of reproducing some results of general rela- 
tivity on the cosmic scale, as well as in deriving the low 
energy Newtonian limit [H, H, 0, HI • The possibility that 
branes can naturally produce a solution to an unsolved 
problem in gravity, such as dark matter, will generate a 
great boost in the theory aside from being a significant 
achievement and a good verification of the branes and 
extra-dimensions ideas. 

In this paper we analyze the weak-field perturbations 
around a flat background generated by matter on the 
brane. We do so in the framework of unified brane grav- 
ity 0; following Dirac's prescription of careful variation 
in the region of the brane (loj (we present the basic prin- 
ciples in the next section). After some remarks on the 
general scenario, we focus on the radial case, thus study- 
ing the field created far from a source. Our main results 
are as follows: 

1. We recover a Newtonian 1/r potential. 



2. The conventional (cosmological) Newton constant 
is suppressed by a constant factor greater than 
1. This difference between cosmological and radial 
Newton constants gives rise to natural cosmological 
dark matter. The amount of dark matter, charac- 
terized by the ratio between the two constants is 
an arbitrary parameter at this point. 

We later discuss the possibility of a transition scale be- 
tween the two Newton constants. Such a transition will 
result in an effective deviation from the Newtonian po- 
tential. An observer who is unaware of this transition 
could interpret it as a continuous distribution of dark 
matter. To calculate the exact transition, one would 
need to analyze weak-field perturbations around a cos- 
mological brane, which is very complicated. We do cal- 
culate roughly the scale of the above transition without 
the exact solution. Remarkably the scale ~ 10 5 ly (light 
years) is only 1 order of magnitude above the experimen- 
tal scale. It remains to be seen whether this prediction 
is correct and whether this transition will lead to fiat 
rotation curves. 



2. The Basics of Unified Brane Gravity 

Dirac has shown [l^] that when performing variation 
of action on a surface, around which one or more of the 
fields are discontinuous, it is crucial to perform the vari- 
ation in a coordinate system where this surface remains 
static, to preserve the linearity of the variation. Vio- 
lating this principle results in nonlinear variation, and 
if this problem is untreated, it would lead to incorrect 
equations of motion. Dirac demonstrated this in his pa- 
per, where he performed both the naive variation and the 
correct variation on a bubble model of an electron, and 
showed that the naive variation results in a missing term 
in the equations of motion. In [ll[, Karasik and David- 
son demonstrated how the naive variation would lead to 
a wrong Snell law, whereas the Dirac prescription leads 
to the correct equation. 

Unified brane gravity^ is based on the same action 
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principle as the standard brane models (Randall- 
Sundrum, Dvali-Gabadadze-Porrati, and Collins- 
Holdom) but following carefully Dirac's prescription for 
correct variation of the brane. We work in a coordinate 
system, where the variation of the bulk metric on the 
surface of the brane has only 5 degrees of freedom, 

Sg ab = 9AB,cSy c y A a y B b + 2g AB y A Jy B b , (1) 

where gAB is the bulk metric and y A are the bulk coor- 
dinates. This way, the brane remains undeformed during 
the variation. The other degrees of freedom are not lost 
but are simply expressed by the constraints that define 
the brane. The equations of motion remain covariant and 
are independent of the reference frame. Using this princi- 
ple of variation, the unified brane gravity field equations 
for a Z2 symmetric AdS bulk with an AdS scale 6 _1 take 
the form 

5 (2) 
3b 1 / 1 \ { ' 

In addition to the familiar terms (namely, the Israel junc- 
tion term, the brane surface tension, the Einstein tensor 
associated with the scalar curvature 7^.4, and the phys- 
ical energy-momentum tensor — SjZ mat ter/Sg >il ' of 
the brane), unified brane gravity introduces A M „. The 
latter consists of Lagrange multipliers associated with 
the fundamental induced metric constraint g tn/ (x) — 
gAsiyix^y^y 3 ,. In the above field equations, X^ serves 
as a geometric (embedding originated) contribution to 
the total energy-momentum tensor of the brane. If the 
variation would be performed naively, it would yield 

A"" = . (3) 

This is reminiscent of Dirac's missing term. A vanishing 
A M „ in Eq. ((2|) results in the familiar Collins- Holdom equa- 
tions (Dvali-Gabadadze-Porrati with AdS bulk), where 
these in turn contain the Randall-Sundrum and Dvali- 
Gabadadze-Porrati equations as special cases. Following 
Dirac's variation principle, we find that X^ is nonzero 
but is conserved, and its contraction with the extrinsic 
curvature vanishes 

A% = , X^K^ = Q. (4) 

This is reduced to the Regge-Tietelboim theory in the 
static bulk limit. Since the Regge-Tietelboim theory was 
derived in a static bulk and is the simplest of all brane 
models, it is very reassuring to have it as a limit (the 
standard brane models are unsuccessful in that). 



3. General Perturbations and the 
Graviton 

We begin with the simplest scenario of a four- 
dimensional flat brane of positive tension embedded in 
five-dimensional AdS bulk 

ds\ = dy 2 + e- 2b ^r lllu dx^dx v . (5) 

b^ 1 = — 6/A5 denotes the AdS scale, r]^ is the four- 
dimensional Minkowski metric, and the brane is conve- 
niently located at y = 0. Before turning to the main 
discussion concerning perturbations of this brane, it is 
imperative to understand the full potential of the unper- 
turbed brane. In the conventional Randall-Sundrum and 
Collins-Holdom scenarios, in order to ensure its flatness, 
the brane has to be of positive (or negative) tension: 



Unified brane gravity, although it requires the same, al- 
lows for one more degree of freedom. 

To see the point, first recall the unified brane gravity 
field equations (|2I4|) . For a flat brane embedded in a five- 
dimensional AdS background, which is the special case 
of interest, K^ v = —brj^. In turn, Eq.(|3]) simply im- 
plies that the corresponding A M „ is traceless. A traceless 
and conserved source serves as an effective (positive or 
negative) radiation term. 

The flatness of the unperturbed brane can be achieved 
the conventional way, if the energy-momentum and the 
embedding terms both vanish, that is, = X^ v = 0. 
But now there exists the milder option 7),„ + X^ = 0. 
Following the above, if (and only if) the real matter 
on the brane exclusively consists of radiation, one can 
choose an appropriate X^ v to cancel it out. To be more 
specific, let our unperturbed flat brane host a constant 
radiation density p, and choose the embedding countert- 
erm to be X% = -T°<™ d = -diag (p, \p, \p, \p). This 
reflects the peculiarity that a flat brane can in fact be 
hot, which is unique to unified brane gravity. The per- 
turbations around such a brane are expected to be quite 
different from those around a Collins-Holdom brane, thus 
giving rise to new physics. To study the perturbations 
induced by an arbitrary source 8T^ V = t^, we find it 
useful to invoke Gaussian normal coordinates, such that 
SgAB — h^S^Sg are the only allowed nonzero compo- 
nents (and reserve the option of supplementing this gauge 
later by the traceless nontransverse gauge). It is impor- 
tant to keep in mind that, although is arbitrary, the 
perturbations of the metric are accompanied by built-in 
perturbations of all of the brane components that are per- 
turbed not by the source itself but rather by the shift in 
the brane space-time structure. For example, the radia- 
tion energy-momentum term must still satisfy the conser- 
vation and traceless conditions, but it must be satisfied 



3 



in the new metric. To this extent, the radiation term is 
corrected via a perturbation T™ d = T°„ rad + 6T™ d that 
satisfies 



d u 5T rad = 



lt\ /"per rpQ^i 



liv 1 crrirad i[ivrpO,rad 

'I " J )iv — n ± iiv J 



ptr rpO,rad\ 
l \v 1 l j,rx ) 



+ h vX d x Tl; ad 



(7) 

(8) 



to preserve conservation and tracelessness, respectively, 
l 

2 



Here T* u = \ (-d x h^ + d^hl + d u h*) is the affine con- 



nection. The above perturbation does not represent an 
addition of radiation (which can be present indepen- 
dently via t^ v ) but rather a geometric effect. By the 
same token, the embedding term is also perturbed via 
A M „ = X° + SX^l, and satisfies 



ft is " 



(9) 
(10) 



However, since for a general perturbation SK^ is not 
proportional to W , the term 
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is not necessarily zero. One can furthermore verify that 
is conserved and not necessarily traceless: 



= n^"? — — A 
— '/ ffiv — 25 
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+ 2b 



(12) 



The nonlocalized part of the perturbation equations is 
the same as the familiar Randall-Sundrum case, since 
the bulk still follows the normal five-dimensional Einstein 
equations 
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d\y\ 



(13) 



where IT] = rf^d^dy is the four-dimensional (unper- 
turbed) d'Alembertian. The localized part of the equa- 
tion is 
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(14) 



The propagation of modes into the bulk remains the same 
as in all of the familiar cases. Thus, we will be focusing 
on only the perturbations on the brane. Performing sep- 
aration of variables, /i M „ = A{y)h lul (x^) 12], where we 
have normalized without loss of generality A(0) = 1 and 
A'(0) 

define a = 1 H — — . Next let us separate the perturba- 



tion hnu = hpfi + h)?J to the standard term K^u 1 , which 
follows the usual brane equation and thus admits the fa- 
miliar solutions and the new term h$ , which is a direct 



2b 



» 
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result of the additional effective source s MJ , . For h 
can write 
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we 
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(15) 



where Grs — bG$ is the Randall-Sundrum gravitational 
constant on the brane, whereas for the new term 



ab 2 



(16) 



Unfortunately, we cannot find a general Green function 
to Eq. p^|) . because there is no closed- form expression 

(u) 

of s^ v in terms of tifaj . To that end, the only general 
prescription to solve Ea.lfTB]) is perturbatively in p (see 
Appendix A). Despite not being able to find a general 
solution, we can get a clue on its properties by taking 
the trace of Eq. (|14[) and reorganizing the various terms 
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(17) 



(11) Keeping in mind that X^ u ~ — p and Gn ~ bG$, by 



looking at the first term in the equation, one may expect 
that the effective Newton constant may take the following 
form: 
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where 
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G 



us 



P_ 

b 2 ' 
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G4 



(18) 



(19) 



is the effective Newton constant in the Collins-Holdom 
scenario and (3 is a dimensionless constant. In the next 
section, we show that this prediction is indeed true and, 
interestingly, (3 is geometry dependent. 

Although we did not obtain a propagator for the gravi- 
ton, the form of the equation looks all too similar to 
the usual brane equations and along with Eq. lfTTj) sug- 
gests that, despite deviating from the standard Randall- 
Sundrum scenario, the graviton propagator remains the 
same. In the following section, we show that the New- 
tonian potential is recovered for large r and thus prove 
that the graviton's zero-mode is massless. 

4. Static Radial Source 

In all studies of gravitational perturbations, the point- 
like radial source is of special interest. Since an exact 
radial solution is missing in all brane theories, the best 
idea we have for a radial potential comes from pertur- 
bative treatment. We solve the equations far from the 



source, in the region where t„ 



0. 
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For the radial case, we show that an exact (non per- 
turbative in p) weak-field solution can be obtained. This 
is mainly due to the fact that we are able to express s MI/ 
explicitly. We choose to work in a traceless Gaussian 
frame. For a radially symmetric perturbation, we can 
choose, in addition to the Gaussian traceless gauge, the 
radial gauge [l3| . It follows that sgg — s^ip = 0. Solving 
the conservation equation for along with Eq. (|12|) and 
gauging following the above, we have 



stt(r) = s rr (r) = ~s(r) + ^ 



drrs(r). (20) 



From Eq. (|16p . we see that since Su = s rr it follows 

This is the familiar form of 



that h 
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radial fluctuations. Finally, hffi constitutes the famil- 
iar Collins- Holdom solution. The exact solution is quite 
complicated, but to first order in 1/r, the solution simply 
yields 

IGghM 



/j( m ) 



(21) 



where M = J d 3 xT U is the mass of the source. Substi- 
tuting Eq. (|12l2QI21j) along with 
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(22) 



into Eq.([Tn]), we can write the equation for hS u > 
K\rhW" + 4/^/ l («)" 

+ {k- 2 -a P yY {U) ' (23) 



-2kh (u) 



AGcnMap 



where Ka 



and k = 



3r 
ab 2 



The solution of 
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physical relevance is the nonhomogeneous one, namely, 



2G CH M 1 
3b 2 r 



(24) 
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The full perturbation h^ u = hffi + hftj is therefore 



htt — h rr — 
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A-kGrsp 
3b 2 



(25) 



It is important to note that it is only due to the solu- 
tion being independent of a that we can proceed without 
integrating over all the mass modes. The Newtonian po- 
tential is thus recovered, giving us further reassurance 
that the graviton is indeed massless, since a mass term 
in the propagator would have generated an exponential 
decay. The associated Newton constant is 

Gch 



Gr 
JV 



1 



AttGrsP 
3b 2 



(26) 



where the r index stands for radial. 

Now that the mathematics has been understood, we 
return to physics. Alone, Eq. (f26|) has nothing new to 
offer. However, gravitational measurements in our Uni- 
verse, although they began with the Solar System, which 
is physically a radially symmetric system, are now quite 
based in the field of cosmology as well. We recall (see Ap- 
pendix B) the cosmological result for expansion around 
a flat background gives an FRW solution with an associ- 
ated Newton constant 
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Grs 



1 

G?4 
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(27) 



where the c index stands for cosmological and p here has 
the exact same role of background radiation. Equation 
(f26l) can also be written as 
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(28) 



Now, if we further assume that the role of radiation in 
our case is also played by the background radiation from 
cosmology, we can compare the two results. First of all, 
since we do have bounds on b from both particle and 
gravitational localization, we can clearly state that the 

term — r is negligible in both equations. This means that 
b z 

G% 



1 N 



Gch, whereas 
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(29) 



The last term in the radial gravitational constant would 



have been negligible if not for the factor 



Grs 

Ga 



We have 



no experimental or theoretical bounds on the latter ratio. 
In fact, the proposed self-accelerated Dvali-Gabadadze- 
Porrati solution for the cosmological constant requires 
this quantity to be very large. If it is large enough, then 
the above term can be significant in the calculation of 
the Newton constant. Thus, in principle, we have a real 
difference between the cosmological and the radial grav- 
itational constants, the radial constant being necessarily 
lower. However, historically, the Newton constant was 
measured in radial systems (the Solar System). Thus an 
observer that is unfamiliar with this physics would in- 
terpret this effective growth of the gravitational constant 
as missing cosmological mass (since, in general relativity, 
mass is inseparable from the gravitational constant), thus 
bringing him to the phenomenon of cosmological dark 
matter, without facing dark matter in the Solar System. 

Although we have not shown it here (this is a conjec- 
ture subject to future research), when solving the pertur- 
bation equations around a cosmological background, one 
expects the two branches of the solution, one being the 
G r N and the other G%, to be connected, creating some 
sort of transition between them. Such a transition, to 
an observer that is unaware of this effect, will seem as 



5 



a gradual increase of mass, that may result in flat rota- 
tion curves. Although the exact solution to fluctuations 
around a cosmological brane is highly complex, we can 
give a rough estimate to the typical scale of such a tran- 
sition. We assume the scale to be roughly in the region 
where the cosmological and radial curvatures are of the 
same order of magnitude, so that the cosmological and 
radial solutions "mix". The radial curvature is of the 
order — , r s being the Schwarzschild radius and the cos- 
mological curvature is of the order of H 2 , H being the 
Hubble constant. The scale of the predicted transition is 
therefore 



r,in 



(rv/ 



.1/3 



Hubble) 



X 



(30) 



where tnubbie is the age of the Universe. When this scale 
is calculated for the Sun, the result is 100 ly, which is way 
beyond the scale of the Solar System. At these distances, 
other stars contribute, and thus the effect is unmeasur- 
able today. For a galactic mass, on the other hand, the 
result is of the order of 10 5 ly, which is only 1 order of 
magnitude higher than the real galactic scale. One needs 
to remember that it is only a rough estimate and also 
that galaxies are not radial systems and are composed of 
many stars, each giving an effect on the scale of about 
100 ly, so that the combined effect may be closer than 
the above result, to give the exact scale of flat rotation 
curves. 



5. Summary and Conclusions 

We have studied the behavior of weak-field perturba- 
tions around a flat brane, in the framework of unified 
brane gravity. It was shown that, even for the most gen- 
eral perturbation, the novel embedding term is "harm- 
less" and the graviton propagator is intact, leaving the 
graviton massless. We verify this result, in particular, for 
a spherically symmetric source, where the conventional 
Newtonian potential 1/r is recovered. However, upon 
a closer examination, we see that, although the func- 
tional form of the potential is standard, the gravitational 
constant differs from the one found in cosmology. Fur- 
thermore, the radial gravitational constant is necessarily 
lower than the cosmological one. For an observer, famil- 
iar only with Einstein's general relativity, this would be 
immediately interpreted as cosmological dark matter. 
This can also be the source of galactic dark matter. The 
flat rotation curves may simply represent the transition 
between the radial and cosmological gravitational con- 
stants. The scale of the suggested flat rotation curves is 

1/3 

predicted in this case to be of the order of (r s t 2 age ) 
When evaluated for a galactic mass, this is indeed close 
to the galactic scale. Despite this transition being the 
natural outcome of the two different gravitational con- 
stants, there is no reason why such a transition would 



generate flat (rather than some general form) rotation 

curves, and the flatness of the rotation curves is wishful 

thinking at this point. 

Although the radial dark matter solution is completely 

speculative in this paper, the cosmological dark matter is 

fully postulated. The only thing that is arbitrary is the 

amount of dark matter. This is due to the arbitrariness of 
Grs 



Gd 



In fact, in order to account for the right amount 



of dark matter, we would need an extremely large G5, 
implying a very low five-dimensional plank mass M5 « 



Appendix A: Perturbative Method 

We can expand the solution to Eq. p^|) via 

00 

( 31 ) 

i=i 

and 

(sab +5350)^-18. m 



where, for i > 1, 



and 



*W = S W^ = aX^hV-V (33) 



S W = 8$rT = a^hM. (34) 



Appendix B: Cosmological Gravitational 
Constant 

In Q we have proven that, when expanding the cosmo- 
logical equations around a flat background with positive 
tension and radiation density of Eq. (76) , 



pip) 



2 lo 



(35) 



where w is a constant. The resulting FRW equation was 
given by Eq.(81): 
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where e = 3- — % — and, therefore, 
er 
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We would like to express the last term in Eq. ([37)1 in terms 
of p and b and, therefore, 
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